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Abstract 



^ ' We provide a fully general-relativistic treatment of cosmological perturbations in a universe 

, permeated by a large-scale primordial magnetic field using the EUis-Bruni gauge-invariant for- 

' malism. The exact non-linear equations for general relativistic magnetohydrodynamic evolution 

, are derived. A number of applications are made: the behaviour of small perturbations to 

' Friedmann universes are studied; a comparison is made with earlier Newtonian treatments of 

, cosmological perturbations and some effects of inflationary expansion are examined. 

O^. 1 Introduction 
u 

bX)' The study of cosmological magnetic fields has many motivations. Any primordial magnetic field 
could provide a seed field for dynamo amplification in disc galaxies, and fields of order 10"^'' gauss 

K/i \ would create significant magnetic fields in galaxy clusters through adiabatic compression alone, [||- 

^ ■ @]- It would also introduce new ingredients into the standard, but necessarily uncertain, picture 



of the very early universe. Large-scale magnetic fields introduce anisotropies into the expansion 
dynamics. They would not survive an epoch of infiation although it is conceivable that large-scale 
fields and magnetic inhomogeneities with a constant curvature spectrum could be generated at 
the end of infiation. A number of proposals of this sort have been made, |9|-|l^]. They involve 
speculative changes to the nature of the electromagnetic interactions in the universe at the time of 
infiation which, as yet, have no other motivation other than to allow magnetic field generation to 
take place. While such changes are not altogether unreasonable, and may reveal important aspects 
of fundamental physics, it is hard to test them independently. Ways of generating magnetic fields 
at cosmological phase transitions have also been explored by a number of authors, [p^]-|16|. 



Any primordial magnetic field must be consistent with a number of astrophysical constraints 
upon its strength in the early universe. Since it provides an additional form of relativistic energy 
density at the epoch of cosmological nucleosynthesis, it increases the expansion rate of the universe; 
hence the neutron-proton freeze-out of weak interactions occurs at higher temperature, with a higher 
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value, leading to an increase in the synthesised abundance of helium-4, [|l7|-p0|]. Therefore helium-4 
observations (extrapolated to zero metallicity) lead to an upper limit on the energy density of any 
cosmological magnetic field at the epoch of nucleosynthesis. If the field is spatially homogeneous 
over large scales then the COBE microwave background data provide the strongest limits on the 
magnetic field strength at last scattering of the microwave background photons, ||2l|, |^^. The 
strength of this limit arises from the subtle evolution of cosmological magnetic fields in expanding 
universes during the radiation era. There is a non-linear coupling between the magnetic density and 
the anisotropic distortion of the expansion needed to support it. The anisotropic pressure created 
by the magnetic field dominates the evolution of the shear anisotropy and it decays far more 
slowly than if the pressure is isotropic, |23|, |22|. This shear distortion determines the microwave 



background anisotropy directly. If the field is inhomogeneous, then the limits weaken, and on small 
scales the fiuctuations can be dissipated. A number of discussions of the damping processes have 
recently appeared, [24|, |2^], and there is a possibility that a field of sufficient strength might have 
observable effects on the pattern of Doppler peaks expected in the microwave background on small 
scales, [^], but the field strengths necessary may not be compatible with the nucleosynthesis and 
COBE limits |2|. 

In this paper we shall focus upon the problem of establishing a formalism for studying the 
general-relativistic evolution of magnetic inhomogeneities in an expanding universe. We do this 
by extending the gauge-invariant treatment of cosmological perturbations, introduced by Ellis and 
Bruni for perfect fiuid cosmologies, to the case of electric and magnetic fields in a universe which 
also contains a perfect fluid. This formalism identifles a combination of gauge-invariant variables 
which specify the evolution of the perturbations in an invariant manner. It offers several advantages 
over the gauge-invariant approach of Bardeen, most notably by virtue of its transparent physical 
interpretation. 

The study of small perturbations to an homogeneous and isotropic Friedmann-Robertson- 
Walker (FRW) universe is beset by subtleties of interpretation because of gauge dependencies. 
The problems arise because two different spacetimes are employed: a 'smooth' spacetime, W, cor- 
responding to the unperturbed background universe, and a 'lumpy' spacetime, W, that is 'close' to 
W both dynamically and kinematically, and represents the 'perturbed' universe we live in. These 
two spacetimes must be related in a way that permits the unique recovery of W from W and vice- 
versa. This is impossible when the only information provided is that the 'distance' between them 
is small in some suitable sense p^ . Instead, one must establish a one-to-one map, ^ : W ^ W, 
so that every point, P, in W has a unique image, P, in W and vice-versa. Establishing such a 
correspondence, that is, a 'gauge' (or a diffeomorphism) between the two manifolds W and W, is 
equivalent to introducing a fictitious smooth background model into the universe. The perturbation 
in the density fi at some point P oiW will be = fJ^ — JI, with Jl the density at P = $~^(P) in W. 
Clearly, there exists a relation between the gauge choice and the value of the perturbation itself. 
Any variation in the correspondence between W and W, which keeps the background model fixed, 
is called a gauge transformation, and it is to be distinguished from a coordinate transformation 
which merely relabels events in the two spacetimes. A gauge transformation not only induces a 
coordinate transformation, but also changes the point in W that corresponds to a given point in 
W. As a result, the perturbation value will depend on the gauge choice if the perturbed quantity is 
non-zero and its spacetime coordinates will depend on the background universe [^]. Even scalars 
which depend on spacetime coordinates, such as the density itself (i.e. /x = /x(t) in W), are affected 
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by changes in the gauge. This creates a problem, since by varying ^ one could give the perturbation 
any desired value at any particular point of W [^] . 

There exist natural choices of gauge (e.g. fundamental fluid flow lines in both W and W, hyper- 
surfaces of simultaneity and homogeneity in W). Unfortunately, some residual freedom remains, 
which is enough to create the aforementioned arbitrariness in the perturbation value. We could deal 
with this gauge freedom simply by keeping track of its consequences, but this is rather inefficient. 
Alternatively, we could specify the gauge completely at the beginning and perform all calculations 
in the chosen gauge. In practice however, difficulties emerge because different gauges are best suited 
for different applications. A third approach, which side-steps the gauge problem altogether, is to 
tackle cosmological perturbations by means of quantities that are entirely independent of the way 
the two spacetimes W and W are associated with each other. 

A covariant gauge-invariant analysis of density fluctuations was introduced by Ellis and Bruni, 
and initially applied to a dust-dominated almost-Friedmann universe |27|. Their geometrical ap- 
proach is simple and physically transparent. In contrast to the standard approach, which compares 
the evolution of the perturbed quantities along the worldlines of two observers, one in W and the 
other in W, brought together via some specific gauge, the covariant gauge-invariant technique com- 
pares the evolutions of neighbouring observers in the same universe, W. It provides exact, fully 
non-linear, propagation equations for all variables, which can be linearized about a "variety" of 
background models. This formalism has already been extended to the study of density irregulari- 
ties within almost-FRW spacetimes with pressure [^], [|^], or with a multi-component fluid [31|, 
as well as to an almost-Bianchi type I universe |3^. Here, we shall develop the covariant gauge- 
invariant technique for a universe that contains a perfect fluid and large-scale magnetic and electric 
fields. This provides, for the first time, a fully general relativistic treatment of electromagnetic 
perturbations for use in cosmological investigations. 

The outline of our analysis is as follows: in sections 2 and 3 the covariant formalism is intro- 
duced. In section 4, the key gauge-invariant variables are defined. Three of them, representing 
spatial changes in the energy density, the pressure, and the expansion within an almost FRW 
universe, were first introduced in ||2^. The fourth gauge-invariant variable is new: it determines 
variations in the magnetic field between neighbouring fundamental observers in almost Friedman- 
nian space-times. The medium, a single perfect fluid of infinite conductivity, is specified in section 
5, together with the relevant equations that characterize its behaviour. In section 6, Einstein's 
field equations are employed to derive the fully non-linear exact propagation formulae for the basic 
kinematic quantities and for all the gauge-invariant variables. Some aspects of the spatial geometry 
are briefly discussed in section 7, and in section 8 we linearize the general equations about a FRW 
universe containing a homogeneous magnetic field. Such an assumption is a valid approximation 
provided that the field is "weak", with an energy density much smaller than that of the isotropic 
perfect fluid, which does not destroy the model's isotropy. Our results are compared with earlier 



non-relativistic treatments by Ruzmaikina and Ruzmaikin |33] and by Wasserman When the 
pressure terms are removed from our equations, the standard results of the Newtonian approach 
appear naturally. We consider several particular applications of the formalism: the case of a dust 
universe; the effect of a magnetic field on a period of de Sitter inflation; and the growth of isocur- 
vature perturbations. 
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2 Spacetime Splitting 



Assume that there exists a well defined average velocity vector in the universe. As a result, space- 
time is supplied with a preferred vector field Ui corresponding to a congruence of worldlines, known 
as the fundamental fluid- flow lines, that carry a family of special observers, namely the fundamental 
observers.^ As usual, the velocity 4- vector is normalized so that UjU* = — c^. 

The existence of the timelike vector field Ui, defined by the motion of the matter in the universe, 
generates a unique splitting of spacetime. To every observer there corresponds, at any instant, a 
tangent spacelike three-dimensional hypersurface T,±, orthogonal to Uj. Such 3-surfaces are the 
observers' instantaneous rest-spaces and they are generally different from each other. All these 3- 
spaces mesh together to create a single hyperplane, the common rest-space for all the fundamental 
observers, only in a non-rotating universe. 

The metric of is provided by a second-order symmetric spacelike tensor, namely the ob- 
server's "projection tensor", defined by 

hij = Qij + -^UiUj, (1) 

and satisfying 



V = K'hi, (2) 



K = 3. (3) 



The projection tensor, together with the velocity 4- vector, decompose tensor fields and tensor 
equations into their spatial and temporal parts. They split the covariant derivative of Ui into 



irreducible basic kinematic quantities, |^], |36] 



0, 1 ■ /.N 

Ui-j = aij + LOij + —hij 2"*""^' 

o C 

where aij = hfj^h.'^u^.^q — @hij/3, cuij = h^^h.^u^.^q, Q = n*.j and Ui = Ui-ju^ are respectively the 
shear tensor, the vorticity tensor, the expansion scalar and the 4-acceleration. The first three terms 
in comprise the second order tensor Vij = h^^h^Uk-^q, which in turn determines the relative 
spatial velocity between neighbouring worldlines. The expansion scalar is used to introduce a 
representative length scale {S) along the observer's worldline. In particular, we define 

for the volume expansion. 



^Hereafter, an observer will always be a fundamental one unless stated otherwise 



4 



The Weyl tensor, Cijkq, is decomposed into a pair of symmetric traceless and completely space- 
like second-order tensors, Eij and Hij, respectively known as the "electric" and the "magnetic" 
parts of Cijkq 13, m, 



Eij = \Cikjqu''u'^, (6) 



c 

H^J = ^V.''Ckq,sU^u\ (7) 

where, rjijkq is the covariant permutation tensor of the spacetime. 

3 The Electromagnetic Field 

3.1 The Electromagnetic Field Tensor 

The electromagnetic field is represented by the second-order antisymmetric Maxwell tensor Fij, 
satisfying 

E[ij-k] = Eij-k + Ejk-i + Fki-j = 0. (8) 
The components of Fij are derived from a potential, Vi, according to 

where Vi = (V^, —^/c) are the covariant components of the 4-potential relative to a freely- falling 
frame. The quantities and ^ are the vector and the scalar potentials respectively. 

3.2 The Electromagnetic Field Components 

As seen by an observer moving with 4- velocity Ui, the electromagnetic field tensor decomposes into 
an electric (Fj) and a magnetic (Hi) part, defined by |36], 



E, = F,ju^ = -Fiju\ (10) 

and 

Hi = Y^V^JkqU'F>^1, (11) 

respectively. From definitions (|To|), ( [ll| ) and the total skewness of Fij and rjij^q, we obtain 
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Eiu' = 0, 



(12) 



Hiu' = 0, (13) 

to ensure that both Ei and Hi are 3-vectors in the observer's rest space. 

The projection of Fij onto the 3-surface the observer's instantaneous rest space, provides 
the decomposition 

Fij = ^ {uiEj - EiUj) - ^r]ijkqu''H'', (14) 

and shows how Ei and Hi completely determine Fij. 

3.3 Maxwell's Equations 

The first of Maxwell's equations has the form 

F% = -J\ (15) 

where Jj is the 4-current which generates the electromagnetic field and obeys the conservation law 
J*.j = 0. The second equation, 

F[^j■k] = 0, (16) 

is a direct consequence of the existence of the 4-potential. 

As measured by a fundamental observer equipped with the projection tensor hij, equations (p^) 
and (ll^) decompose into temporal and spatial parts as ||3^ 



E\i + 2jHi-\E'ui = ec, (17) 



a) + uj) - ^/i^j E^ + ^r]'^'"^UjUkHg - crf^^'iujHk.q - ch'^J^ = h'jE^, (18) 

H\i-^jEi-^Wui = 0, (19) 

+ Jj - ^/i^) W - ^n'^^iu.UkE, + -J^'^'^u.jEk,, = h'^H\ (20) 

where the vorticity vector (tOi), and the charge density (e = —J^Ui/c^) are new variables. Note that 
LOij and LOi are related by the formula 
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^rjijkgUJ^u'', (21) 



which ensures that LOjUJ^^ = 0. 



3.4 The Energy-Momentum Tensor of the Electromagnetic Field 

The electromagnetic energy- momentum tensor (T*^) is a symmetric and trace- free second-order 
tensor of the following general form 

= F''F^ - \g'^F,,F'^. (22) 



Relative to a fundamental observers, T*^ splits as , 

^el = ^ + + ^ (5 + + ^u^'v'^'"^'ukE,H, + M^^, (23) 

where E'^ = EiE^ and = HiW are the magnitudes of the field's electric and the magnetic 
components respectively; Mij is a traceless and completely spacelike symmetric tensor given by 

M'^ = -(^ + hA h'^ - \e'E^ - WW. (24) 

The expression (^^ allows a direct comparison to be made between TJ^ and the energy- 
momentum tensor, *Tij, of a general imperfect fluid possessing viscous and heat conduction con- 
tributions: 

= yu'u^ + *ph'^ + ^*q^'u^^ + *7r'^. (25) 



The resulting correspondence, which provides a "fluid description" for the electromagnetic field, is 

Mem = ^(5 + ^')' (26) 

I ( e"^ \ 

Pem = -l^-^ + H^y (27) 

qL = U'^'u.EkH,, (28) 

<in = (29) 

and leads to the familiar equation of state of radiation = ^^c? /?>. The last expression suggests 
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that Mij contributes an anisotropic electromagnetic pressure. 



4 The Key Gauge-Invariant Variables 
4.1 The Di, Yi and Zi Spatial Gradients 

In their covariant and gauge-invariant study of cosmological perturbations Elhs and Bruni 
describe spatial variations in the energy density (/ic^), the pressure (p) and the expansion (0) 
by projecting their gradients onto the instantaneous rest space of an observer comoving with the 
expanding fluid. Assuming that the unperturbed background universe is represented by a FRW 
spacetime, they consider the following basic variables 



(30) 
(31) 
(32) 



where k = SirG/c^ is the Einstein gravitational constant. All three vectors Xi, Yi and Zi vanish 
in a FRW model, thus ensuring their gauge-invar iancej^ Ellis and Bruni then define the following 
"comoving" variables 

A = (33) 

namely the comoving fractional orthogonal spatial gradient of the energy density, and 

Zi = SZi, (34) 

which is the comoving orthogonal spatial gradient of the expansion. Besides being invariant under 
gauge transformations, Di and Zi also describe the spatial variations of //c^ and G respectively, 
within a perturbed FRW universe 



4.2 A Gauge-Invariant Variable for the Magnetic Field 

The assumption of a large-scale primordial magnetic field requires that the fictitious background 
universe must include some degree of anisotropy and a spatially homogeneous non-rotating Bianchi- 
I spacetime provides the simplest example of such a universe. However, later we shall primarily 

^The simplest gauge-invariant quantities are scalars that are constant in the background, or tensors which are 
zero there. The only other possibility is a tensor written as a linear combination of products of Kronecker deltas with 
constant coefficients (Stewart and Walker lemma [S]). 
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consider the case of an unperturbed FRW spacetime. This is an acceptable simphfication provided 
that the field is too weak to affect the model's isotropy and allows us compare our approach with 
the existing (Newtonian) results directly. 

Although Ellis and Bruni restricted their study to almost-FRW geometries, the gauge-invariance 
of their key variables also holds in any spatially homogeneous, non-rotating background spacetime, 
such as a Bianchi-I universe. Furthermore, the same variables can describe spatial gradients in fic^, 
p and Q within an almost Bianchi-I model, as it has been argued by Dunsby [02| . 

Let us consider a spatially homogeneous universe permeated by a large-scale magnetic field (i.e. 



Hi = Hi{t)) with no accompanying electric field (i.e. Ei = Qi). Following Ellis and Bruni |27], we 
define the second-order tensor 

n^, = Kh,^h^Hk,q = K(3)Vj/7i, (35) 

namely the orthogonal spatial gradient of the magnetic field. When the vorticity is zero, Tiij van- 
ishes (see Appendix A.l), and so provides a new gauge-invariant variable. Moreover, the quantity 

Mij = SHij, (36) 

namely the comoving orthogonal spatial gradient of the magnetic field, describes the spatial vari- 
ations of Hi between neighbouring world lines in a perturbed FRW universe (see Appendix A. 2). 
Following Dunsby [ |32| , we can extent this result to a nearly Bianchi-I spacetime if required. 
The quantity TCij is a completely spacelike and traceless tensor 

u'Hij = HijU^ = 0, (37) 
n\ = 0, (38) 

where the last result is derived from (|l9|) under the assumption of a zero electric field. 

5 Specifying the Medium 

5.1 The Case of Infinite Conductivity 

In order to specify the material content of the universe, we consider the covariant form of Ohm's 
law HI, 

Ji + \j^UjUi = -aEi, (39) 

where Ui is the fluid velocity and a represents the conductivity of the medium. Projecting onto the 
instantaneous rest space of a fundamental observer, we obtain 
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hi J 3 = Z^^i- (40) 



1 

— ( 

c 

Thus, non-zero spatial current densities (i.e. /i/ Jj 7^ 0), are compatible with a vanishing electric 
field as long as the conductivity of the medium is infinite (i.e. a 00). Under this assumption, 
formulae ([l7D-(pO[) become 

2jHi = ec, (41) 
i]'^''%[ukHg-c^Hk.g) = c^hy^, (42) 
H\-\wui = 0, (43) 



= h^^W. (44) 



Equation ( |4^ ) can be rearranged as 

h^^H\j = ^^^ViW = 0, (45) 

to provide the familiar 'vanishing 3-divergence' of the magnetic field. Some useful expressions can 
be obtained from (^); in particular, 

/i^ (^S^W) ' = (ct^ + w^) S^W , (46) 

is the covariant analogue of the 'induction' equation and verifies that within an exactly FRW 
universe the magnetic field decays adiabatically as the inverse square of the scale factor. An 
additional formula, which will be used later to simplify the energy-density conservation law, is 
obtained after contracting (Q) with Hi 

The introduction of a perfectly conducting medium, as opposed to the assumption of a pure 
magnetic field with no electric field and zero spatial currents, is important in this approach. As 
we shall see, it allows perturbations in the energy density of the fiuid to be coupled with magnetic 
irregularities in a straightforward and natural way (see further discussion in Appendix B). 

5.2 Magnetized Perfect Fluid with Infinite Conductivity 

In a medium with infinite conductivity (i.e. Ei = 0), the electromagnetic energy-momentum tensor 
(see (|23|)) takes the form 
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TiL = ^u'u^ + ^h'^ + M'^, (48) 

where 

M'^ = —W^ -WW. (49) 

The energy-momentum tensor of a single perfect fluid (i.e. zero energy-flux and no anisotropic 
stresses) is 

T}^ = fiu'u^ +ph'^. (50) 

Thus, the energy-momentum tensor that describes a single magnetized perfect fluid of infinite 
conductivity has the form 

T'l = l^/i + ^ j u'u^ +(^p+^^ h'l + M'^, (51) 

with trace T = 3p — fic^. 

A comparison between (51) and (p5|), the stress tensor of an imperfect fluid, shows that a single 
magnetized perfect fluid of infinite conductivity can be represented as an imperfect fluid with the 
following properties 



W 

V = /^ + ^' (52) 
b 

*q = 0, (54) 

*TT^^ = M'^ = —h'^ - WW. (55) 

We can now derive the conservation laws that characterize the medium. We start by setting the 
divergence of (|5l|) equal to zero. The resulting formula is then split into a temporal and a spatial 
part, respectively expressing the conservation of the energy and of the momentum densities. More 
precisely, contracting with the observer's 4-velocity and using (^) we find 

fic^ + (^c^ + G = 0, (56) 

which coincides with the energy-density conservation law for a single non-magnetized perfect fluid. 
The conservation of the momentum-density is obtained after projecting ( [5l| ) onto the observer's 
instantaneous rest space. 
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and can be rewritten, in order to involve the spatial gradients Yi and fiij, as 

f^ + ^ + ^jUi + Y,- 2n[ij]W + ^u^Mji = 0. (58) 

Equation (^) (or (|58|)) naturally connects spatial fluctuations in the magnetic field, with the 
acceleration and gradients in the pressure, and subsequently with spatial variations in the energy 
density (see also Appendix B). 



6 The General Propagation Equations 

6.1 The Field Equations 

The general form of Einstein's field equations is 

kT 

J^ij = ^Tij _ _g^l + ^glj ^ (59) 

where Rij = Rj^jj^ is the Ricci tensor, gij is the spacetime metric and A is the cosmological constant. 
The Ricci scalar (R = ii*-) and the trace (T) of the energy-momentum tensor, that describes the 
material content, are related by 

R = AA- kT. (60) 
Using the observer's projection tensor, the Ricci tensor splits as 

Rij = h\h\R^i - ^ {h\R'i^UkU^ + u'ukR^'ih^) + ^R'^^UkUgu'u^ , (61) 

which shows that the decomposition is entirely determined by the sums h^j^h^^R'"^, h^jR^^u^ and 
R^^UiUj. 

When dealing with a single perfect fluid of infinite conductivity, permeated by a large-scale 
magnetic field, Tij is given by ( [5l| ) and consequently we find 

h'kh'.R"'^ = (^|(^/ic2-p+:|^^+A^/^*^■ + KM^^ (62) 
h'jR^'^Uk = 0, (63) 

2 

R^^UiUj = ^ (;uc2 + 3p + H^) - kc^. (64) 
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Thus, the Ricci tensor associated with our cosmological model, decomposes into purely spatial 
and temporal parts only, 

R'^ = (f (/.c^ - ^ + f^) + ^) + + (i + I + ^ ) - I) (65) 
while the Ricci scalar becomes 

R = {^iic' - 3p) + 4A. (66) 

These equations allow us to derive the formulae that determine the time-evolution of the ex- 
pansion, the shear and the vorticity. 



6.2 Propagation of the Kinematic Quantities 

The evolution of basic kinematic variables such as the expansion scalar (0), the shear tensor {(Jij) 
and the vorticity vector (i^j), along the world line of a fundamental observer, is governed by their 
respective propagation equations. These are derived from the following formula, describing the 
evolution of the relative velocity between neighbouring fundamental fluid-flow lines 

hl'h^Vkq - ^UiUj - ^'^'^^jUi + VikV^j + RikjqU^vf = 0. (67) 

Contracting (|67[), recalling that Vij = Uij + tOij + @hij/3, and using (|^), we find that under the 
influence of a cosmological magnetic field Raychaudhuri's relation becomes 

e-A + 2[a'-u;') + ^ + ^ (/.c^ + 3p + H^) - Ac' = 0, (68) 

where A = ii^.^, a' = and to' = by definition. 

The shear propagation formula is obtained by taking the symmetric, trace- free, part of (|67|). In 
order to derive the final expression, we employ (^), the definition of the Weyl tensor along with 
equations (|62[), (64) and (|6^). The outcome is 



2R .. 

— d*^ + w^w^ + = 0, (69) 
3 

where Eij is the electric part of the Weyl tensor and Mjj is given by (^9|). In deriving (69) we 



have used the result Wj^w^- = WjWj — io'hij, which relates uj' to the vorticity tensor and vector. 
Comparing (|69| ) with the corresponding formula in the case of a non-magnetized imperfect fluid 



(see [36]), we notice again that the distortion generated by the field is analogous to that induced 
by anisotropic stresses. The proportional to the projection tensor term, simply subtracts off the 
trace. 
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Finally, the skew part of (^) , provides the equation governing the rate of change of the vorticity 
vector 

h'^ [S^J) ■ = ^Ji^'iu.uk-, + (70) 

where, = rf^^'^UjOJkq/'^c. Notice that the magnetic effects are felt indirectly through the shear 
and the acceleration. 



6.3 Propagation Equations for the Key Variables 



From definitions (^)-(|3J) we can derive the propagation formulae for the various spatial gradients. 
The most important equation determines the evolution of the comoving fractional orthogonal spatial 
energy-density gradient (Dj) 

-^n. + -,um,^. (71) 

It is obtained by taking the time derivative of (|3^) and then using the conservation laws (|56| ) and 
(|5^. Notice that the field gradients are already connected with those in the energy density, unlike 
the Newtonian approach, where the coupling occurs at the next level. This coupling arises from the 
presence of the quantity iii, which will then bring Di and M^j together via the momentum-density 
conservation law. 

The vector Zi, the comoving orthogonal spatial gradient of the expansion, has a rate of change 
determined by the following propagation formula. 



3kS 

T 



STZiii + ^iimji + SAi - 25(^)Vi [a 



2 Co' 



(72) 



where Ai = h/ A.j = (^^VjA, and 



3c^ c' 



UJ 



A 3 



(73) 



with the quantity in brackets representing, as we shall see in section 7, the 3-Ricci scalar of the 
observer's instantaneous rest space. To derive (72) we start with the time derivative of (34), and 
then use (|58|), together with (38). 
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The third key variable is the comoving orthogonal spatial gradient of the magnetic field {Mij). 
In order to derive the propagation formula for Mij, we first split the covariant derivative of the 
field with respect to irreducible kinematic quantities, employing the observer's projection tensor 
and equation to obtain, 



1 2 2 1 

Hi-j = -Hij ^Hka'^^-Uj] + —HkU)^^jUj) + {2HiUj + UiHj) -vf'HkUiUj. (74) 

The time evolution of the spatial gradient of the magnetic field can now be obtained by applying 



the Ricci identity to the field vector, in connection with (74). The outcome, 

S-^h,''hf (s^Mk,) ■ = -M,k {a) + io)) + [a,'' + Mk, - '^H^Z, - 

nS ( B \ 

-^u^Hk + + -KA - KSh^^RkqjsH'^u^ (75) 

shows how the Riemann curvature tensor Rijkq acts as an additional source for magnetic inhomo- 
geneities. Notice that so far there has been no formula determining the rate of change of spatial 
gradients in the pressure. When an equation of state for the non-magnetic fluid is introduced, the 
propagation of Yi, along the observer's worldline, can be obtained from ([7l]). 



Dunsby |31[] applied the covariant and gauge-invariant technique to the study of cosmological 
perturbations in an almost FRW universe filled with an imperfect non-magnetized fluid. As men- 
tioned earlier, our stress can be written as corresponding to an imperfect fluid with the special 
properties given by (|5^)-(|55|). When Dunsby's fully non-linear propagation formulae for Di and Zi 
are applied to this particular imperfect fluid, we recover equations ( [7l| ) and ([7^). Of course, there 
is no expression equivalent to ( |75|) in Dunsby's analysis. 



7 Aspects of the Spatial Geometry 

The three-dimensional Riemann tensor {^^^ Rijkq), that determines the curvature of a fundamental 
observer's instantaneous rest space, is defined via the commutator for the 3-D gradients of any 



spacelike vector. In particular, provided that fjU* = 0, we have |3C], 



(3)v,(3)v,-t;fe - (=^)v/=^)Vi^;fe = -^i0^jh^'Vq + ^''^ Rqkj^V^ (76) 



where, by definition 



Rijkq = hfhfhf^h^Rpsrt - \vikVjq + \viqVjk, (77) 
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with Vij = ^^^VjUi = aij + Uij + G/ijj/3. The successive contractions of ( ^7|) provide the 3-Ricci 
tensor i'-^^Rij) and the 3-Ricci scalar of E_l, by 

^^^Rij = hi'^h-'Rkq + ^Rikjgu''ui + - ^Vij, (78) 



<»)B = if + |ft,uV-^ + |!-^. (79) 



and 



Notice the vorticity term that appears in equation (|7q). In a rotating spacetime, this term prevents 
the commutation between the 3-gradients of scalars. In particular, the formula 

(3)vX3)v^./ - (3)V/3)V./ = (80) 



holds for any scalar quantity /. This non-commutativity, reflects the fact that the fluid velocity 
does not consist a hypersurface orthogonal vector field if ujij ^ (see Appendix in [^] for more 
details) . 

Equations ( [78| ) and (^) characterize the curvature of the observer's rest space within a general 
spacetime. We shall rewrite ( |78| ) and (|7^ ) for a magnetic universe containing a single perfect fluid 
of infinite conductivity. Using (|6^), ( [67| ) and (^), we obtain 



Moreover, from (p4\j, ( pq ) and ([79|), we find the following 3-curvature scalar 

K.I»)ii = 2(.(M.nf:)-g + ^-^ + A). (82) 

Combining ( |8l[ ) with (p2|), we can now decompose the 3-Ricci tensor into its trace and its trace-free 
parts 



\h,''hf(S-''(S'u:kJ -U[k;g]). (83) 



The evolution of /C, along the observer's worldline, is governed by the following propagation 
equation, derived from (H), {Mj, (|6|) and iM) 



16 



(/C - I (a^ - ■ = (/C + + f (a' - u^) - 2.a.,WK (84) 

The 3-curvature scalar qualifies as an additional gauge invariant variable, provided that it is con- 
stant in the unperturbed cosmological model. In a spatially homogeneous irrotational background, 
the 3-gradient /Cj is always independent of the gauge choice, as is implied by the relation 

. ... ^ (a + ^,M,H^) - 5^,2. + I'" V. (.^ - . (85) 

which will be useful in discussing the case of isocurvature perturbations. 

It should be emphasized that all the relations given above refer to the instantaneous rest space 
of an individual fundamental observer. For non-rotating spacetimes, the same formulae (without 
the vorticity terms) determine the geometry of the observers' common rest space. 



8 Approximations 

8.1 Linearization About a FRW universe 

Let us assume that the background universe is described by a FRW spacetime containing a perfect 
fluid and a weak large-scale magnetic field (i.e. <^ /^c^). The gauge-invariant variables are 
the shear (cjjj), the vorticity (coij), the acceleration (ttj), the divergence of the acceleration {A), its 
spatial gradient {Ai), the electric and magnetic parts of the Weyl tensor {Eij, Hij), the spatial gra- 
dients of the energy density (Xj, Di), the spatial gradient of the pressure (Yi), the spatial gradient 
of the expansion (Zj, Zi) and of the magnetic field {Ti.ij, Mij), together with the spatial gradient 
(/Cj) of the 3-curvature scalar. All these quantities vanish in a Friedmann universe. According to 
(p^, the 3-curvature scalar is independent of the gauge choice, provided that it vanishes in the 
background spacetime. 

The formulae above are the exact, fully non-linear, propagation equations. We shall linearize 
them about the FRW model. In the process we shall retain the mass density (//), the pressure 
(p), the expansion (G) and the magnetic field {Hi), which are all spatially independent in the 
background, as zero-order quantities. Variables that vanish in the FRW background, together 
with their derivatives, will be treated as first order. The anisotropic pressure {Mij) generated 
by the primordial field will be also considered as of order one, since it must disappear within a 
FRW spacetime. Finally, terms of order higher than the first will be disregarded. Under these 
conditions, the conservation law ( |56|) for the energy-density remains unaffected, whereas law (^), 
for the momentum-density, changes. In particular, we have 

fic^ + (jic^ + e = 0, (86) 

and 

K(^x + ^Ui + Y,- 2W[ij] W =Q. (87) 
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Raychaudhuri's equation, (|68|), becomes 



e - ^ + ^ + ^ (^fic^ + 3p) - Ac^ = 0, 



3 2 

and the time-evolution formula for the vorticity vector, ([70|), reduces to 

Linearization permits the omission of the projection tensor from the left-hand side of the remain- 
ing propagation equations and the weakness of the field implies that Di/2 + MjiH^ / k^c? ~ Di/2. 
Hence, equations (71), ( |7^ ) and (|75|), which describe the evolution of Dj, Zi and M.ij respectively, 
become 



(90) 



2G Kuc^ 
= -—Zi - -h^Di - 3c^M[ij]W + STZiii + SAi, 



(91) 



where now 7^ = /C/2, with /C = 2 (k/zc^ - eVSc^ + A), and 
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2k 



-^u^Hkhij - nSh.^RkqjsH'^u'' . 
Linearization leaves Maxwell's equations unaffected: hence, 

- \wui = 0, 



{2HiUj + iiiHj) + 



(92) 



(93) 



and 
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(94) 



The key role played by the acceleration (uj), by its divergence {A = u^.j) and by the latter's 
spatial gradient {Ai = ^^^Vj^) will become clear below. Therefore, it is useful to derive expressions 
for these three variables. The momentum-density conservation law (^7|), immediately provides a 
useful formula for Uj, that leads to the following expression for A 



A 



(95) 
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This is derived by applying the commutator (|76D to the spacehke vector Hi, and substituting the 
zero-order part of the Ricci tensor associated with from (|83|). Furthermore, the spatial derivative 
of (H), combined with (||), (H), (^) and (|8|) provides an expression for Ai within a perturbed 
FRW universe: 




/C,-2ef4 + :77^^^1^'^V,^/' (96) 



3(// + p/c^) Ic^ 3(^c^+p) 



where = p/ [jl is the sound speed and Bi = ^^^Vji^^ is an additional first-order variable. 

We make one further comment about the 3-curvature scalar /C. Its linearized propagation 
formula, 



^ = -f(^ + |^), (97) 



contains a first-order term generated by the divergence of the acceleration. Thus, to zero-order, we 
find 

= (98) 



where k is the 3-curvature constant (i.e. A; = 0) of the background FRW spacetime [27|. One can 



replace /C by (|98|) whenever the 3- Ricci scalar is coupled to a quantity of first or higher order. Such 
a substitution further simplifies equations ( |9l| ) and (96). We shall adopt this approach in what 



follows. Finally, equation (85) reduces to 



2K^l^ 49 
= ~~S~ ~ 35c2 ( 

which can be used to modify (p6|), if required. 
8.2 Pressureless Fluid 

Let us assume that the perturbed FRW universe, described in the last section, is dominated by a 
single non-relativistic (i.e. p = 0) perfect fluid. Under this restriction the equation of continuity, 
), becomes 

fi + ^e = 0, (100) 



and leads to the familiar evolution formula for the mass density of a dust-dominated cosmological 
model 
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_ ^ 



(101) 



with M = 0. 

An important result emerges from the conservation of the momentum-density. In particular, 
the contraction of the pressure-free form of (^) with the magnetic-field vector provides the relation 

UiW = 0. (102) 



Hence, within a magnetized dust-dominated nearly FRW universe, the acceleration must always 
be normal to the magnetic field, a result consistent with the Lorentz force law of special relativity. 
Furthermore, equation (|102| ) may be rearranged to give 

Hiu' = 0, (103) 



and ensures that the field's derivative actually lies in the observer's rest space. Also, p = in (p7|) 
leads to the following expression for the acceleration 

U^ = —n[ij]W. (104) 

Kjjjj 



This demonstrates how the geodesic motion of the matter is disturbed by the spatial variations 
of the magnetic field. We can express (|104D in a more familiar form: Ti-^ij] is an antisymmetric 
3-tensor, the generalized curl of Hi projected onto the observer's rest space. Such a tensor has only 
three independent components and corresponds to a 3-vector, namely the covariant spatial curl of 
Hi. In particular (see Appendix C.l), we find 

^fo] = -'^eijkCurlH\ (105) 



where eijk = ijijkqu'^/c is the covariant spatial permutation tensor. When substituted into (104), 
(|105| ) provides the following expression for the acceleration, familiar from classical magnetohydro- 
dynamics 

u^ = -e,,k{{curlW)H^). (106) 

Notice that results ( 102]) , ( |103D and ( 104 ), or equivalently (106), also hold within a pressure-free 
almost Bianchi-I universe. 

An additional familiar result is obtained when we adapt the vorticity propagation formula, (^) , 
to our cosmological model. From a;* = —curlu^/2 and equation (|106D, we find that 



{curlu'y + '^curlu' = -e'^^^'^^V j {ekgs{curlH'^)H') , (107) 
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in agreement with Wasserman's analysis |34]. 

The adaptation of Maxwell's equations to this cosmological model gives 



= 0, 



(108) 



and 



—h\, \w =w 

3 J 



(109) 



with the latter providing the standard expression H'^ oc S*"^ for the evolution of the energy density 
of the magnetic field with the expansion scale factor S{t). 

We shall now proceed with some additional simplifications. Suppose that spacetime is flat (i.e. 
Rijkq = 0) and that there is no cosmological constant (i.e. A = 0), then Raychaudhuri's equation 
(p^) reduces to 



@ = A + 



(110) 



since /C = 2 (k/^c^ — 0^/3c^) .Using (|102|) and (|104|) , the propagation equations for the spatial 
gradients (pO|), (|9l|) and (p^) become 



= — 



26 



M 



(111) 
(112) 
(113) 



respectively, where the pressure- free expression for the quantity j4j, that appears in (|112D , is 



3//c' 



,2 ^3^1 ' 



(114) 



Bi = (3)Vii?2^ and Ki is given by (||). 

The antisymmetry of means that the indices i and j in (111) must take different values 

{i ^ j). Thus, irregularities in the energy density are not influenced by the component of the field 
that acts parallel to them, a result which is also apparent in the analysis given by Ruzmaikina and 
Ruzmaikin [33|. Combining equations ( |lll| )-( p!l4 ) and ignoring the non-linear terms, we obtain 
the following second-order differential equation for the time-evolution of the comoving fractional 
orthogonal spatial gradient of the energy density. 
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" 2B ■ Kuc^ ^ S / 6k |"^^„9^ „ 6k ^ , 

In a perturbed Einstein de Sitter universe (i.e. A; = and A = 0) this reduces to 

26- Kfic'^ 5.o^ 2 465^2, 3^ 2G5.o^ ■ • , , 



The first three terms on the right hand side of (|llq ) also appear in the Newtonian treatments 
presented by Ruzmaikina and Ruzmaikin [33|, and by Wasserman |^^; the first two are immediately 



identifiable, while the third is obtained after a transformation (see Appendix C.2). The latter 
becomes obvious, when the gauge-invariant approach is applied within the non-relativistic context. 



as it has been shown by Ellis in |40]. The resulting second-order differential equation contains 
only the three aforementioned standard terms. Consequently, our relativistic approach provides 
the classical equation plus two extra relativistic terms. The first term links the vorticity to the time 
evolution of linear density gradients through the 3-divergence of the vorticity tensor. It appears 
because in a general spacetime 3-D surfaces need not be orthogonal to the fluid-flow. However, 
the term's coefficient suggests a weak effect since <C /ic^, or no effect at all if the universe is 
static. Such a 'correction term' does not appear in the Newtonian treatment, because there the 



hyperplanes orthogonal to the fluid flow are always tangent to the 3-surfaces of absolute time [30|. 
Also, using (|2l|) and the covariant definition of a vector's spatial curl we find 

495^2, 495/72 
3/ic2 ^ ' 3fic^ ' ^ ^ 

where = —curlu^/2. 

Let us focus on the last term in ( |116| ). It emerges from the coupling between Di and Mij in ([7l|). 
We see that its presence also depends on the universal expansion, but it also creates only a weak 
effect, since it contains the ratio H"^ / ^c? <C 1. Similarly to the acceleration case above, the quantity 
^^^V^jHij represents the generalized curl of the field's derivative projected onto the instantaneous 
rest space of a fundamental observer. Thus, it can be mapped onto the covariant spatial curl of 
Hi. As a result, the second relativistic 'correction term' can be reshaped into the following easily 
identifiable form 

29 •? 9 •? 

Thus, according to (117) and ( |118| ), equation ( |116| ) becomes 



A = -^A + ^ A + ^^'^V^B, + ^^curlu:, - '^e.^,{curlW)HK (119) 
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8.3 Long- Wavelength Solutions 

The gradient Di describes spatial variations of the energy density orthogonal to the fluid flow. One 
can extract the information contained in Di by considering the following local decomposition [30|, 



S^^)ViD^ = A,j = + + ^Ahij, (120) 

where = A(jj) — Ahij/3, Wij = A[jj] and A = S^^^V^Di are all gauge-invariant quantities. The 
symmetric, traceless tensor, Sjj, describes the evolution of anisotropic structures (e.g. pancake 
or cigar-like structures). On the other hand, the skew tensor Wij characterizes the rotational 
behaviour of Di. Hence, the scalar A contains the information regarding the spatial aggregation of 
matter. When dealing with the problem of structure formation, the latter quantity is the crucial 
one [41 1, [42|. The propagation formula for A is found by taking the 3-divergence of ( |115| ) and 



using the properties of the 3-D gradients as well as the fact that, in the linear approximation, the 
total divergence of the vorticity tensor is zero (see section 7 and also [^]). We obtain 

■■ 29: K^c\ fAk ^,^_2^ f 3k 2Q'^\fAk fs^_2^ ..2 / s 



where B = (^^V'-Bj = ^^^V'^H'^. The final form of (|121| ) is obtained only after treating the field's 
energy density (i?^) as a perturbation itself]^ 

Assuming that spatial and temporal dependences are separable, we express every gauge-invariant 
first-order quantity in ( |121| ) as a sum of time-independent scalar harmonics Q^"'\ defined as eigen- 
functions of the Laplace-Beltrami operator. 



(3)v2qW = _^qW, (122) 

where the eigenvalue n directly corresponds to physical wavelengths only if A; = [^]. Thus, in an 
almost Einstein de Sitter universe, and for low frequency (i.e. n — > 0) fiuctuations, we may ignore 
all the terms in the right-hand side of ( |12lD . Consequently, the large-scale evolution of matter 
perturbations is determined by 

AH + - = 0, (123) 

3 2 

where A^") is the nth harmonic component of A, with ('^^VjA^"') ~ 0. Its solution, 

^Although it has not been explicitly used yet, the characterization of as a gauge-invariant perturbative term 
1 line 

/ is very small. 
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= aV"V2/3 + aL"V-\ (124) 

contains a growing and decaying modes, exactly as in the case of a dust-dominated non-magnetized 
universe and in agreement with ||3^ . The variable r measures the observer's proper time, while the 

(n) (n) 

quantities A_|^ and A_ remain constant along its worldline. 



9 Two Special Cases 

9.1 The False Vacuum Assumption 

The introduction of a medium that obeys the "inflationary" false vacuum condition (i.e. p + = 
0), implies that the spatial gradients Xi and 1^, defined by (|30|)-(|3l|), are related by the simple 
formula 



Yi = -Xi. (125) 

The exact propagation equation for the orthogonal spatial gradient of the energy density, Xj, is 
obtained by taking the time derivative of (]30|). From the conservation laws (|5^ ) and (|5^), we have 



5-^ V (s^Xj) ■ = -X^ (a{ + - K [fic^ +p)Zi- 2en[,j]W + 



3,2 + (126) 

where Zi has been defined in (|3^ . In an almost-FRW magnetized universe, the false vacuum 
assumption, together with the momentum-density conservation law, reduces (126) to the following 
evolution formula 

= |, (127) 

with X = 0. Hence, spatial gradients in the energy density die away as S~^, independent of their 
wavelength and in line with the cosmic "no-hair" theorems for expanding inflationary universes. 

According to (]5^), however, the energy density of the material component remains constant 
along the observers worldline. Consequently, definition (^), along with ( |127| ), suggests that the 
fractional density gradients do not change on comoving scales 

Di = const. (128) 



In other words, fundamental observers experience no variation in Di within their event horizons: 
density contrasts remain "frozen-in" as long as the equation of state p = —^(? holds. 
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9.2 Isocurvature Perturbations 

We define isocurvature perturbations as fluctuations evolving in a non-rotating (i.e. uJij = 0) 
spacetime of constant spatial curvature (i.e. /Cj = 0). Within an almost-FRW universe, equation 
) provides a condition for the occurrence of isocurvature inhomogeneities: 
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X, = ^Zi. (129) 

The consistency condition for this kind of disturbance is obtained by linearizing the time deriva- 
tive of ( |129| ) about the background FRW spacetime, giving 

6/c / e \ 2G , 

where Ai is given by (p^). 

Provided that the smooth model has A; = 0, we can ignore the first term in ( |130D , which means 
that isocurvature disturbances cannot survive in a non-static universe so long as Ai ^ 0. On the 
other hand, in the absence of pressure, formula (|9^ ) implies that the remaining term is unimpor- 
tant on large scales. Consequently, the matter era seems capable of preserving long-wavelength 
isocurvature fluctuations. 



(130) 



10 Discussion 

We have pursued the study of cosmological density perturbations in a universe that contains a large- 
scale primordial magnetic fleld by means of the Ellis-Bruni covariant and gauge-invariant method. 
We have deflned a new variable that is independent of the gauge choice and describes the spatial 
variations of the magnetic fleld. Assuming a universe that contains a single perfect fluid of very 
high conductivity, we have derived the exact, and fully non- linear, equations that determine the 
model's time evolution. When the linearized equations are applied to the simple case of pressure- 
free matter, familiar results are recovered. The energy density of the matter is found to evolve 
as the inverse cube of the scale factor and the energy density of the magnetic fleld as S~'^. The 
acceleration is provided by the well known expression of classical magnetohydrodynamics, and is 
always normal to the magnetic fleld in accord with the Lorentz force law. 

We also identify the general-relativistic corrections to the Newtonian treatment. These are 



introduced by the two extra terms (117) and ( 118 ) on the right-hand side of the propagation 



equation, (116), of the density gradient. The flrst term is due to universal rotation and does not 



affect the gravitational clumping of matter. A similar vorticity term also appears in the treatment 



of a perfect fluid with non- vanishing pressure, |3C]. Here the fluid pressure is zero, but there exists 



a residual isotropic pressure, induced into the model by the magnetic fleld (see (p^)). Technically, 
the second relativistic correction, which modifles the evolution of both Di and A, appears because 
everything is projected onto the observer's instantaneous rest space (where all measurements take 
place). Physically, its source is the anisotropic pressure, also introduced by the fleld (see (|55|)). 



In Dunsby's imperfect fluid analysis, |31], such a term is incorporated into the time derivative 

((3)V,7r/);fcn^ 
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By focussing upon the spatial aggregation of matter, we have obtained long-wavelength so- 
lutions for (|115D , which reveal the relative unimportance of the magnetic field on the growth of 
large-scale density fluctuations. We also provide solutions during a period of de Sitter inflation. 
The results show that density gradients decay, but slower than in the non-magnetized case inves- 
tigated by Ellis and Bruni, . Furthermore, within the immediate neighbourhood of a comoving 
observer perturbations freeze out. Our analysis permits a transparent approach to isocurvature 
inhomogeneities. We have obtained a criterion for the occurrence of this kind of disturbances as 
well as a consistency condition. The latter shows that a dust-dominated almost-FRW universe, can 
sustain long- wavelength isocurvature perturbations. 

A major advantage of the Ellis and Bruni technique is that the non- linear propagation equations 
can be often linearized about a variety of background models. For example, we can linearize our 
exact equations about a smooth Bianchi-I universe. This cosmological model is anisotropic. The 
shear no longer behaves as a first-order variable and the model's evolution becomes more compli- 



cated, 1 22]. The formalism we have established here, will enable a full gauge-invariant analysis to 



be carried out to determine the effects of cosmological magnetic and electric fields on the early 
universe, and on the temperature anisotropy of the microwave background. 
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APPENDICES 

A The Spatial Gradient of the Magnetic Field 
A.l Proof of the Gauge-Invariance for J^ij 

The metric associated with an exact FRW or with an exact Bianchi-I spacetime can always be 
written in the diagonal form gij = diag{gQQ, gn, g22, 933) with the quantities ga (no sum) functions 
of time only (i.e. gtj^^ = 0).f\ As a result the associated spatial Christoffel symbols vanish 

= 0. (131) 

Following definition (^), we split the spatial gradient of the magnetic field as 

Khth^H^.fl + Khi%!'Ho.g. (132) 
^Latin indices run from to 3, whereas Greek ones from 1 to 3. 
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In a spatially homogeneous universe iJj^^ = 0. Also, in a comoving frame, equation (|l^ ) means 
that Hq = 0. Also, when there is no rotation, Uj = —c5^, so = 0. Thus, via ( |131| ), relative to a 
comoving frame, equation ( |132| ) gives 



7i 



0. 



(133) 



Therefore, definition (35) provides a spatial gradient for the magnetic field that satisfies the re- 
quirements for gauge invariance, when the background universe is a FRW or a Bianchi-I spacetime. 



A. 2 Description of Magnetic Spatial Variations by M.ij 

Let us consider two neighbouring fundamental observers P and P'; {x*} is a coordinate system 
comoving with the expanding fluid. Relative to this frame, the timelike worldlines associated with 
the two observers, which are assumed to be close enough to measure the same proper time (r), are 
labelled by (x'^ = const, ct) and by (x* + 5x^,ct) respectively. With respect to a general frame 
yi — y^(^x^)^ the vector (5x^^,0), which connects the two worldlines at all times, becomes 

5v^ = ^6xr (134) 

As far as the observer P is concerned, the relative position vector between the two events is the 
projection of Sy'^ onto its own 3-D instantaneous rest space T,± 

S±y' = h'jSyK (135) 

In order to examine the fleld's spatial variation, the observer at P must parallel transport its 
own vector {Hi)p to P' and compare it to the one defined there by the field itself. The resulting 
difference should then be projected onto T,±. The field vector at P' as seen from P is 

{Hi)p, = Hi {y^ + 6^y^) = {Hi)p + {Hij)p 6^yK (136) 

The parallel transport of {Hi)p to P', along the infinitesimal displacement generates the 

vector 

{Hl)p, = {Hi)p + {T%)^{Hk)p 5^yK (137) 

The spacelike part of the difference between ( |136D and ( |137| ) gives the variation of the magnetic 
field at P' as measured by the observer at P. Since P is arbitrary, we have 
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(138) 



Within an almost FRW spacetime, the relative position vector evolves as 5±y'^ 
{^±y^)Q = const, [27|. Hence, we arrive at the formula 



Si Hi 



1 



SH 



{S±y' 



S (5±y*)o! where 



(139) 



and conclude that the second-order tensor Mij = STiij determines the spatial variations of the 
magnetic field between two neighbouring fundamental observers. 



B The Case of a Pure Magnetic Field 

The adoption of an infinitely conducting medium has led to the omission of the electric field from 
Maxwell's equations, while preserving the spatial currents. An alternative approach is to assume a 
pure magnetic field with no electric field and zero currents. This assumption reduces (|l8|) to 

rj'^'"^Uj [ukH, - c'^Hk-,) = 0, (140) 



and subsequently to 



C 



Superficially, (|141| ) seems to provide a relation between the acceleration and the field gradients. Yet, 
when inserted into ([5^), it causes the magnetic terms to cancel out, thus effectively disconnecting 
inhomogeneities in the field from those in the energy density of the medium. Some coupling between 
Di and Mij can still be achieved via the third term in the right-hand side of (|72|). However, this 
coupling is only significant when ~ /j,c^, but is negligible in the case of a weak magnetic field 
discussed here. 

C Aspects of the Pressureless Case 
C.l An Expression for 7i[ij] 

From the definition of the covariant spatial curl of a 3- vector (see ) and the total antisymmetry 
of eijk we have 
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curlW = e'^'^^^'^SJjHk 



-e''%m- (142) 



[jk]- 

Thus, contracting with the e-tensor, and using the identity e^^ ^psk — 

we obtain 

= -'^^ijkCurlH^ . (143) 



C.2 Comparison with the Newtonian Approach 

In Newtonian theory, and in comoving coordinates, the evolution of the density contrast 5 = dfi/fi 
is determined by, |]3^ , 

By contrast, the covariant gauge- invariant approach provides the fohowing propagation formula for 
the comoving fractional orthogonal spatial energy-density gradient, 

Relative to a comoving frame D = dD/dt, so the correspondence between the first three 



terms in (144) and ( |145D is straightforward, due to the established correspondence of Di and 5 (see 
Appendix A in [^^). Also, recalling that Bi = ^^^Vj-ff^, the fourth term on the right-hand side of 
dBl) is 

in comoving coordinates. Since the tensorial quantity Mij/k describes the infinitesimal spatial 
variations of the magnetic field (see appendix A. 2), the direct correspondence between the fourth 
terms in (|144| ) and ( |145| ) is established. 
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